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A model of disordered spin-Peierls system is considered, where domain walls are randomly dis- 
tributed as a telegraph noise. For this realization of the disorder in an XX spin chain, we calculate 
exactly the density of states as well as several thermodynamic quantities. The resulting physical 
behavior should be qualitatively unchanged even for an XXZ chain, up to the isotropic XXX point. 
For weak disorder, besides a high energy regime where the behavior of a pure spin-Peierls system is 
recovered, there is a cross-over to a low energy regime with singular thermodynamic properties and 
enhanced antiferromagnetic fluctuations. These regimes are analyzed with the help of exact results, 
and the relevant energy scales determined. We discuss the possible relevance of such a disorder 
realization to the doped inorganic spin-Peierls compound CuGeC>3. 



One dimensional quantum spin systems in the presence 
of randomness show unusual and intriguing properties 
(see e.g. Ref. and references therein). For instance, 
it has been shown that the ground state of the 

Heisenberg antiferromagnet with random exchange con- 
stants can be interpreted as a random singlet state, where 
pairs of spins are coupled into singlets with an energy gap 
to the triplet configuration which is weaker for widely 
separated pairs. The uniform and staggered magnetic 
susceptibilities, x an d Xs, have a (Griffith's like) singu- 
lar behavior at low temperature, x ~ Xs ~ l/(T\n T). 
Interestingly, in spite of the singlet nature of the ground 
state, the spin-spin correlation functions are still long- 
ranged. In fact, (S'(r)S l (0)) ~ l/r 2 , for large r (i — 
x,y,z) ffl. These properties are not modified by spin 
anisotropy if, on average, J z < J x = J y . Notice that spin 
anisotropy does not manifest itself in the spin-spin corre- 
lation function with different power law behavior of i = z 
with respect to i = x, y, contrary to the case in the ab- 
sence of disorder. The behavior of the random XXZ chain 
is however unstable towards a finite average dimerization, 
i.e. a finite average difference between the exchange con- 
stants of the even bonds and of the odd bonds. This case 
was recently analyzed by Hyman et al. Q by means of a 
real space renormalization group approach. For a finite 
average dimerization <fi, they find that the spin-spin cor- 
relation functions decay exponentially with a correlation 
length £ ~ \(/>\~ 2 , but the Griffith singularities remain, 
even if weaker. In particular, singularities of the uniform 
susceptibility x ~ T"^ 1 , and the specific heat C v ~ T a , 
where a cx \<p\, are found to persist [0. 

The study of the role of disorder in a spin-Peierls sys- 
tem may be useful to understand the behavior upon dop- 
ing of the inorganic spin-Peierls compound CuGeOa. The 
pure compound is known to undergo a structural transi- 
tion at 14K M, below which the Cu02 chains dimerize 



and a spin-gap opens. However, upon substitution of few 
percent of Cu with magnetic (Ni Q) or non magnetic 
(Zn Jt]||) impurities (as well as replacing Ge with Si 0), 
besides the structural transition, which still occurs close 
to 14K, an antiferromagnetically ordered phase appears 
below a lower temperature Tjy ~ 4K. Moreover, the es- 
timated magnetic moment with 4% of Zn is as high as 
0.2/zs ||. This behavior is quite puzzling. First of all, 
heuristically, one would expect a Neel temperature ex- 
ponentially small in the ratio of the average distance be- 
tween the impurities to the spin-Peierls correlation length 
X S p. At 4% doping, this would imply T N /T S p ~ 0.04, 
inconsistent with the experiment. In addition, one would 
also expect a magnetic moment of the order of the doping 
concentration, not almost an order of magnitude larger, 
as seen experimentally. 

In this Letter, we study a particular realization of a 
disordered spin-Peierls system which does show a large 
enhancement of antiferromagnetic fluctuations, coexist- 
ing on a lower energy scale with an underlying dimer- 
ization. Moreover, this model permits an exact calcula- 
tion of physical quantities for a wide range of tempera- 
ture/energy. 

The Hamiltonian of each chain in the absence of im- 
purities is 

H = £(1 + M-IY) (S*S* +1 + SfSf +1 + AS*S* +1 ) , 

i 

(i) 

where <po is the strength of the dimerization. We assume 
that one impurity releases one spin- 1/2 solitonic excita- 
tion, connecting regions of different dimerization parity 
]l0| . The role of the interchain coupling is to provide a 
confining potential to the soliton, which will be trapped 
within some distance from the impurity p!l| . Moreover, 
the weak link connecting the impurity nearest neighbors 
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(which would be for instance generated by a next-nearest 
neighbor exchange) is approximated to be equal to the 
weak bonds in (|l|). Therefore, the effective Hamiltonian, 
defined now on a chain of one site less, remains the same 
apart from the presence of a domain wall. For a finite 



of this function, the density of states of the fermionic 
model is 



number 



1 1 nip 



of randomly distributed impurities, the ef- 



fective model will therefore be assumed to consist of a 
chain with rn mv sites less, described by the same Hamilto- 
nian Eq. ([[]), but in the presence of randomly distributed 
domain walls. This amounts to take a site dependent 
4>(i), which takes alternatively two values ±0o, jumping 
from one to the other at the (random) position of the 
antiphase walls. We will show that it is possible to cal- 
culate many physical properties of the soliton band which 
is created by disorder inside the spin-Peierls gap, without 
the precise knowledge of the soliton wave functions. In 
Eq. (lj), A = corresponds to the XX chain, while A = 1 
is the isotropic XXX model. On the basis of the analyses 
of Refs. [^PQ, we expect that the behavior at < A < 1 
should be similar to that at A = 0, therefore we will only 
study the latter case, which is much simpler. We believe 
that this approximation gives qualitatively good results 
for all the range < A < 1, especially in view of our 
particular choice of the disorder. By means of a Jordan- 
Wigner transformation, the model can be mapped onto 
a model of disordered spinless fermions. By linearizing 
the spectrum around the Fermi energy, introducing the 
right and left moving components of the fermion field, 
and then taking the continuum limit, the diagonaliza- 
tion of the Hamiltonian amounts to solve the following 
coupled differential equations: 
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d 

'T^XfleM + <j)(x)xLe(x) + lh s XLe(x) = eXRe(x), 

d 

where XR(L)e(x) is the eigenfunction of energy e on the 
right(left) moving field, and we have also considered for 
later convenience a uniform staggered magnetic field h s in 
the z-direction. The dimerization field <j){x) corresponds 
to that introduced in Eq.(|l|), apart from an appropriate 
normalization factor. The equations can be decoupled by 
the following transformation 

U+ e (x) = XRe(x) + iXLe(x), 
U_ e (x) = iXRe(x) ~ XLe(x). 

These two functions are solutions of the Schrcedinger-like 
equations 

(~ + ^{ X ) ± U±e( x ) = Eu± e (x), 

where E = e 2 ~h 2 should be greater than zero. In the fol- 
lowing, we will often use the integrated density of states 
as a function of E, which we will define as N(E) . In terms 



p(e) = 2e 



dN{E) 
~9E~\ t 



(2) 



In the case in which <p(x) is a white noise, these equations 
have been analyzed quite in detail in the context of disor- 
dered one-dimensional Fermi systems |l2|,(n| , or classical 
diffusion of a particle in a random medium (for a review 
see e.g. Ref. EJ]). An interesting anomaly of this prob- 
lem is that, for a zero-average white noise, the E = 
state is extended p^Jl5|JT^ |, and both the localization 
length and the density of states diverge as e — > 0. Quite 
recently, Comtet, Desbois and Monthus |l7j (CDM) spe- 
cialized those equations for a particular disorder, for 
which they have been able to calculate exactly the inte- 
grated density of states N(E) and the localization length 
A(E). Specifically, they assumed a random potential 0(x) 
which takes alternatively two values 0o and 0i at inter- 
vals whose lengths I > are randomly distributed accord- 
ing to the probability densities /o(0 = no exp (— uqI) and 
fi(l) — 7iiexp(— nil) (see also Ref. Q)- This choice of 
4>{x) is particularly suited for studying our problem of 
randomly distributed domain walls. In particular, our 
case corresponds to <pi = ~0o < 0, and no = rii, i.e. to 
an average dimerization = (0oni+0ino)/(no+ni) = 0. 
Nevertheless, we will also discuss the more general situ- 
ation riQ ^ n\, in which case is finite. Moreover, we 
start by taking h s = 0. In the model there are three rele- 
vant length scales, Asp — l/0Oi = l/^o an d h = l/ni- 
As p is the correlation length of the system in the absence 
of disorder, which is the case if for instance lo/h — > oo. 
In this case, the spectrum of the single-particle excita- 
tions (which is symmetric around zero energy) shows a 
gap 20o, and the density of states p(e) has an inverse 
square root singularity at e = ±0o- For generic lo and 
h, the density of states can still be exactly calculated 
within the phase formalism approach fl7j| , and expressed 
in terms of integrals which have to be numerically eval- 
uated. Essentially, the method consists in writing the 
master equation for the joint probability distribution of 
the phase of the wave function and (f>(x), and solving for 
the stationary cc-independent solution. In particular, if Iq 
and l\ are much longer than Asp, i.e. if the gap has the 
time to develop in a region of constant </>(x), the density 
of states still shows a peak at ±<fro, even though states 
are created inside the gap. These states accumulate, in 
a singular manner, as e — > In particular, the density of 
states around zero energy goes like p(e) ~ e 2 ^ 1 , where 
/i = (ni — no)/(2(f>o) is finite. In Fig.l, we draw p(e) 
for e > 0, 0o — 1; an d various no and ri\. For fj, = 0, 
p(e) ~ l/|eln 3 e|. The key feature of our choice for the 
random potential is that, even if the average dimerization 
= 0, i.e. if no = n% — n, the density of states shows 
a pseudo-gap if 0o ^> n (see Fig.l), totally absent for a 
white noise process [181. 



2 



To be more precise, from our numerical results we 
find, similarly to CDM, that the integrated density of 
states N(E) for weak disorder (i.e. both no and m much 
smaller than <j> ) saturates below the pseudogap 0o to a 
value ~ n^ni/ (n® + n\), which is of the order of half 
the average number per unit length of steps of the ran- 
dom potential. The saturation occurs at an energy scale 
which can be identified as the typical effective band- 
width of those midgap excitations. This result physically 
implies that, for weak disorder, the number of states 
generated inside the gap is of the order of the average 
number of domain walls. From the analytical expression 
of N(E), we obtain that ln(E*/(p%) ~ — 2(/>o/(n + m), 
i.e. E* is exponentially small in the inverse of the dis- 
order strength. In addition, it is also possible to cal- 
culate the localization length X(E). In particular, for 
4> 7^ 0, A(0) = l/4>, which implies that the localized wave- 
functions inside the gap have a much longer localization 
length than the spin-Peierls correlation length Xsp- More 
interesting, for = 0, which is relevant for our disorder 
modelization, X(E) ~ | In £7 1 , so that the states close to 
E = are almost delocalized. 

More generally, our model at low temperature/energy 
is equivalent to the models analyzed in Ref. 0] and in 
Ref. [||, for if) — and cf> ^ 0, respectively. The anal- 
ogy can be expected by the following arguments. For 
e < tfio and uq — n\ <C <f>o, the problem reduces to a 
model of weakly coupled spins localized close to each do- 
main wall. As a first approximation, only the exchange 
coupling between two successive spins can be retained, 
which is given by J(r) ~ 0o exp(— r&o), being r the ran- 
dom distance between two domain walls distributed ac- 
cording to nexp(— rn). Thus the model is indeed equiv- 
alent to an Heisenberg chain with randomly distributed 
exchange constants. The probability distribution of J at 
energy scales < 0o can be readily found to be 

and it has to be used as the starting point of the renor- 
malization group flow equations of Ref. (ij . In this way, 
it is possible to recover the same results that we obtain 
by exploiting the exact solvability of our model, thus 
showing not only that the two models are equivalent, but 
also that spin-anisotropy does not really matter jl9| . For 
no 7^ rix, the same analogy works now with the model of 
Ref. Q. More rigorously, the above conjectured equiva- 
lence can be proven by showing that the models have the 
same low temperature thermodynamic properties. 

In our model, we can in fact calculate exactly many 
thermodynamic quantities and find not only the low tem- 
perature but also the intermediate (T ~ <j>a) temperature 
behavior. For instance, the uniform magnetic suscepti- 
bility is given by 



r°° dN 1 

X {T) = 13 dE%— - = — , 

V ' Jo d£2cosh 2 (/V£/2) 

and is plotted in Fig. 2 for the same values of 4>o, no and 
ni as in Fig.l. 

From the asymptotic behavior of N(E) for small E, 
we find that, at low T, %(T) - T 2 ^ 1 , for u ^ 0, 
and - l/(7Tn 2 T) for /x = 0. The latter is exactly 
the result for the random XXZ Heisenberg model. Our 
model thus belongs, at low energy and for // = 0, to 
the same universality class. For all /i's smaller than 1/2, 
the magnetic susceptibility still diverges at low tempera- 
ture. Analogously, the specific heat vanishes as C v ~ T 2p 
(C v ~ 1/| In 3 T\, for /i = 0), which is compatible with the 
result of Ref. [Q with 2/i = a, thus showing the equiva- 
lence with our model at /i ^ 0. In addition, we obtain the 
full behavior of x at intermediate temperatures, as shown 
in Fig. 2. We see that, at T ~ tfio, the susceptibility de- 
creases as if a spin-Peierls gap were present, even though 
it finally diverges at low T. Moreover, for E* < T < <pQ, 
we predict a Curie like behavior, with a Curie constant 
oc AT*. 

The behavior of the staggered part of the spin-spin cor- 
relation function Xs{ x ) can be deduced by the analogies 
with the models analyzed in Refs. Jl],|]]. In particular, 
for (j) ^ 0, Xs{x) decays exponentially with a correlation 
length cx (1/^i) 2 0. On the contrary, for the case rel- 
evant to our model, which corresponds to <f> = 0, Xs{%) 
decays as a power law ~ 1/x 2 0. At finite tempera- 
ture and n — 0, h\x s {x,T) ~ — xa j ln 2 (T/^o) 0, where 
a = 20q/(?t,o + ni). This expression suggests a new en- 
ergy scale E c , which can be identified as the coherence 
energy for the antiferromagnetic fluctuations. In fact, 
when T > E c , the correlation function should behave like 
exp (— 2x(f>o), which leads to ln(£' c /0o) ~ ~V a / ( l ) o, that 
is to a coherence energy exponentially small in the in- 
verse square root of the disorder strength, but still much 
bigger than E*. The appearance of an energy scale gov- 
erning the spin-spin correlation function, which differs 
from that entering the average density of states, is not 
unexpected in the presence of disorder, which introduces 
basic differences between average and typical behaviors. 
On the other hand, for ix ^ 0, below another energy scale 
Efi , we should recover the result of Ref. || , which sets 
ln(£W0 2 ) ~ 

We also exactly calculate the staggered susceptibility 
Xs {T) . By means of Eq. (||) we find that 

*■<?>- f o) 

For u < 1/2, this susceptibility diverges at low T like the 
uniform susceptibility. However, while the integral over 
E in the uniform susceptibility is cut-off by T 2 , the con- 
tribution to the singular behavior of the staggered suscep- 
tibility comes from all E up to approximately E„ . More- 
over, all higher energies also contribute to the staggered 
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susceptibility with a finite term as T — > 0. Therefore, 
while the singular behavior deriving from all e = \f~E~ < T 
can be ascribed to local excitations, that deriving from 
e > T is solely due to longer range antiferromagnetic fluc- 
tuations Q| . The rapid enhancement of antiferromagnetic 
fluctuations that we find is extremely suggestive in the 
light of that recently observed in CuGeC>3, as previously 
discussed. In fact, our model for a disordered spin-Peierls 
system clearly shows a coexistence of dimerization with 
long range antiferromagnetic fluctuations, the latter ex- 
isting on energy scales lower than the pure spin-Peierls 
gap. These fluctuations may induce a magnetic ordering 
via the interchain coupling, below some Neel tempera- 
ture Tm- The magnetic susceptibility would then still 
show the drop at the Peierls transition, but the low tem- 
perature divergence would finally be cut-off by TV, below 
which xCO would exponentially vanish, compatibly with 
the experimental evidences (see e.g. Ref. Moreover, 
on the basis of the previous discussion, we expect that 
the Neel temperature is related to the energy scale gov- 
erning the spin-spin correlation function, that is to E c , 
which is exponentially small in y/l/Xsp, and therefore 
larger than the typical bandwidth of the low energy ex- 
citations E*, which is exponentially small in l/Xsp- This 
difference might be the explanation of the relatively large 
Neel temperatures found in the doped CuGeC>3. 
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FIG. 1. Density of states for (f>o = 1 and no = m = 0.3 (dotted line), no = m = 0.1 (full line), no = 0.1, m = 0.3 (dashed 
line). Also shown in the insert is the low energy behavior. 
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FIG. 2. Uniform magnetic susceptibility at zero staggered magnetic field, for the same cases as in Fig.l 
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